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Abstract 



The surface modes of a composite consisting of aligned metallic wires with 
square cross sections are investigated, on the basis of photonic band struc- 
ture calculations. The effective long-wavelength dielectric response function is 
computed, as a function of the filling fraction. The dependence of the optical 
absorption on the shape of the wires and the polarization of light is discussed, 
and the effect of sharp corners analyzed. The effect of the interaction between 
the wires on the localization of surface plasmons is also addressed. 

Keywords: Interface states; Light scattering; Plasmons; Electron-solid inter- 
actions, scattering, diffraction. 
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The interface electromagnetic modes of composite structures [[!]] are well-known to be 
of basic importance in the understanding of a number of physical problems such as optical 
absorption 0, van der Walls attraction ||, surface-enhanced Raman scattering and the 
response to incident charged particles 0-§[. Analytical descriptions of normal electric modes 
have been limited to slabs and isolated spheres P,|TT| and spheroids |TT|, while numerical 



methods have been required for most shapes |12| including such simple cases as single cubes 
PllUjl and rectangles [0 . 



Dipolar modes in a system of identical spheres embedded in an otherwise homogeneous 
host have been approximately described, for many years, from the knowledge of the effective 
transverse dielectric function first derived by Maxwell- Garnett |15| (MG) within a mean- 
field approximation valid for small values of the volume occupied by the spheres. Also, 
many attempts have been made to account, at large filling fractions, for higher multipole 



interactions, which are absent in the MG theory |[T6|-t20|. Recently, new techniques have 
been developed for solving Maxwell's equations in structured materials pT|-f2^]. Based on 
these numerical solutions of Maxwell's equations, exact calculations of the long-wavelength 



limit of the effective transverse dielectric function of periodic dielectric [21 and metallic 



25^pT| structures have been performed. From these calculations, strengths and positions of 



surface modes of a composite made up of long metallic cylinders have been obtained |27 
and the MG theory has been shown to offer a good approximation as long as the distance 
between the axis of neighboring cylinders is larger than twice the diameter of the cylinders. 



In this paper we extend the work presented in Ref. [^7| to the case of aligned metallic 



wires with square cross sections. We consider a binary metal-dielectric mixture with a 
volume fraction 1 — / of the insulator, with a real positive dielectric constant eo, and / 
of a periodic system of long metallic wires with square cross sections arranged in a square 
array with lattice constant a = xd, as shown in Fig. 1. The cross section of the wires is 
taken to be small in comparison to the wavelength of the electromagnetic excitation and 
large enough that a macroscopic dielectric function e(uj) is ascribable to the wires. Thus, 
the electromagnetic properties of the composite will be defined by a single effective dielectric 



function, 



k 2 c 2 



L0 Z 



(1) 



the wave vector k being that corresponding to a Bloch wave propagating through the com- 
posite. For simplicity, the magnetic permeabilities are assumed to be equal to unity in both 
media. 

The dielectric function of metals is nearly real and negative below the plasma frequency 
down to frequencies comparable to the inverse of the conductivity mean free time. Hence, 
negative values of e(a>) give rise to a variety of surface resonances, which allow to express the 
long-wavelength limit of the effective dielectric function of the composite in the following 



form 28,291: 



e c ff(^) = e 
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where 



u(uj) = [1 - e(u)/e ] 



-i 



(3) 



m v (0 < m u < 1) is the depolarization factor associated with the vth normal mode, and B v 
(B v > 0) represent the corresponding strengths, which all add up to unity: 



(4) 



Similarly |]30 



e fr l^J 



where < n v < 1, C v > 0, and 
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(5) 



Ea = i. 



(6) 



Both depolarization factors and strengths in Eqs. (2) and (5) depend only on the microge- 
ometry of the composite material and not on the dielectric function of the components. 



Optical absorption is directly given by the imaginary part of the effective dielectric 
function of Eq. (2), and coupling to charged particles is described by the imaginary part 
of the effective inverse dielectric function |HJ] of Eq. (5). In a homogeneous metal (/ = 1), 



there is only one mode strength different from zero, with depolarization factors mo = and 
n = 1, where Ree — > ±oo and Ree ~ 0, respectively. Thus, for a metal described by the 
Drude dielectric function, e(u) = 1 — ui^/iv 2 , the bulk resonance at the plasma frequency, 
Up, can only be excited by penetrating charged particles. 

In the case of two-dimensional composite structures (see Fig. 1), homogeneous along 
the axis of the inclusions, there are two different values of e e s(oo) corresponding to different 
polarizations. If the electromagnetic (e.m.) wave incident on the structure is polarized along 
the axis of the inclusions (s-polarization), the presence of the interfaces does not modify the 
electric field and the response of the composite is easily found to be equivalent to that of a 
homogeneous medium having the effective dielectric function of Eqs. (2) and (5) with only 
one mode strength different from zero, B — Cq = 1, and depolarization factors tuq = and 
n = f P2| , |5B| . Thus, for this polarization the energy-loss function, Im[— e^f (u;)], shows a 
single peak at a reduced plasma frequency. For a metal described by the Drude dielectric 
function, the reduced plasma frequency is u = \fjuj p . 

Now we focus on the response to e.m. waves polarized within the plane normal to the 
axis of the inclusions (p polarization). For this polarization all depolarization factors m v 
and n v of Eqs. (2) and (5) are found to satisfy the relation |27| 



n v = 1 - {D - l)m„, (7) 

where D represents the dimensionality of the inclusions, i.e., D = 2. As long as only dipole 
surface modes can be excited, B\ — C\ — 1, a combination of Eqs. (2) and (5) with Eq. (7) 
results, for D = 2, in the depolarization factors mi and rti to be given by 

m 1 = ±(l-f) (8) 

and 



m = ^(i + /). (9) 

Eqs. (2) and (5) with B\ — C\ — 1 and the depolarization factors of Eqs. (8) and (9) 
reduce to the two-dimensional MG effective dielectric function [01] , as discussed in Ref. 
We note that multipolar modes can only be neglected in a few simple cases, in which the 
long wave-length effective dielectric function is accurately given by the MG approximation. 
These are the cases when the two-dimensional structure is made of (a) a sparse (/ — > 0) 
distribution of circular inclusions, and (b) a dense (/ — * 1) distribution of square inclusions. 

For a single inclusion embedded in a host material with dielectric constant eo, both 
Im[e e ff(u;)] and Im[— e^(uj)] are proportional to the particle susceptibility, and satisfy the 
relation 

InfeV)] = Jm[e eS (u))/el (10) 

i.e., the spectral representations of Eqs. (2) and (5) coincide. Thus, one finds from Eq. (7) 
that, in the case of single inclusions, for any given mode with m u < 1/D there exists another 
mode with depolarization factor m' u > 1/D such that 

m' v = 1 - (D - l)m u . (11) 

In particular, for single circular inclusions only dipole surface modes can be excited, depo- 
larization factors being mi — n\ — 1/2, and both the optical absorption and the electron 
energy loss exhibit a single strong maximum at the dipole resonance where u — 1/2, i. e., 
e(a>) + €q = 0. For single inclusions with shapes other than circular, this resonance splits 
into several peaks that occur at a set of frequencies according to Eq. (11). 

In order to give a full quantitative description of all surface modes in our two-dimensional 
composite structure, we have followed the method developed in Ref. |23[ for calculating 



dispersion relationships of Bloch waves in structured materials, and we have computed the 
effective dielectric function from Eq. (1). The insulating component of the composite has 
been chosen to have a real dielectric constant of eo = 1 and the dielectric properties of the 
metal inside the inclusions have been modeled by a Drude dielectric function of the form 



where u p represents the bulk plasma frequency of the metal and 7 is an inverse electron 
relaxation time. The plasma frequency of the conduction electrons in aluminum has been 
used, %uj p = 15.8eV, and the parameter 7 has been chosen to be small enough so as to 
distinguish all multipolar excitations: 7 = O.leV. 

Figs. 2 and 3 show our numerical calculations of Im[e e ff(u;)]// and Im[— e~ s (uj)]/ f, as a 
function of the frequency co, for various values of the ratio x between the lattice constant 
and the side of the wires. For x > 6 our results are almost insensitive to the precise value 
of x, and they nearly coincide with the effective dielectric function of single wires. In this 
(dilute) limit, Im[e e fr(o;)] and Im[— e^u;)] coincide, and they both exhibit several resonances 
satisfying Eq. (11) (see the insets of Figs. 2 and 3) fl35| . 

At higher concentrations of metal, both the optical absorption (Im[e e ff(u;)]) and the elec- 
tron energy loss (Im[— e^f (00)]) show a dominant dipole-like mode and a band of multipole 
resonances. We note (see Figs. 2 and 3) that as x decreases the dipole-like resonance in 
Im[e c ff(u;)]// and Im[— e^jr (a;)]// is shifted towards lower and higher frequencies, respec- 
tively, as a consequence of the interaction between one surface and another. This result, 
already predicted within the MG approximation (see Eqs. (8) and (9)), is also found in the 



case of circular inclusions [27]. Also, it is interesting to notice that as long as x decreases 
multipole resonances tend to be negligible, thus dipolar depolarization factors approaching 
those predicted by Eqs. (8) and (9). 

In the case of circular wires that are touching (/ = 0.785) the metal forms a connected 
medium, and the so-called energy-loss function Im[— e^(o>)] shows, therefore, a bulk plasmon 
excitation at u = uj p (n\ = 1). Besides this peak, there is also a band in Im[— e^j (&>)] from 
multipole contributions to the effective response at uj v = ^Jn u (f)u p [y = 2, n v < 1), which 
prevents our computed effective dielectric function to coincide with the MG approximation 
(see Ref. [^7|). However, the presence of flat surfaces with sharp corners, as in the case of 
square inclusions, results in all multipolar strengths to be negligible for / > 0.75 (x < 1.15); 
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thus, our computed effective dielectric function is well described, in this case and for x < 1.15, 
by the MG approximation. 

Finally, we have represented in Fig. 4 universal curves (solid lines) for the dipolar mode 
depolarization factors and strengths, versus the filling fraction, as determined with use of 
Eqs. (13) and (14) from our computed effective dielectric function of a system of infinitely 
long wires with square cross sections. An inspection of this figure indicates that the trend 
with increasing filling fraction is for the dipolar peaks in Im[e e ff (a>)]/Im[— e~£(u}) to approach 
the MG results of Eqs. (8) and (9), also plotted in this figure (dotted lines); furthermore, 
the MG results nearly coincide with our numerical results for / > 0.8 (x < 1.1). Associated 
with the deviation, for smaller filling fractions, of the actual dipolar mode positions from the 
MG results is the reduction in the dipolar mode strengths Bi and Ci, which are represented 
in the inset of Fig. 4. These mode strengths have been determined by 

B v = 4rf ImkffM] (13) 

and 

C u = ^lm[-e-iM], (14) 

where io v is the frequency associated with the z/th normal mode, and H represents the peak 
height in the bulk energy-loss function [in the case of the Drude dielectric function of Eq. 
(12), m„,n v = ool/ujp and H = lo p /^\. 

In summary, we have presented, on the basis of ab initio numerical solutions of Maxwell's 
equations in structured materials, exact calculations of the effective long-wavelength dielec- 
tric response function of a periodic system of long metallic wires with square cross sections. 
We have investigated the dependence of both the optical absorption and the electron energy 
loss on the shape of the wires and the polarization of light. In the case of s polarization 
the effective dielectric function does not depend on the shape of the wires and there is 
only one mode strength different from zero. In the case of p polarization it had already 



been concluded [[27]] that MG results are good as long as the distance between the axis of 



neighboring circular cylinders is larger than twice the diameter of the cylinders. Now we 
have concluded that in the case of two-dimensional periodic structures consisting of square 
inclusions MG results are reliable as long as the ratio between the lattice constant and the 
side of the inclusions is smaller than x ~ 1.1. Both the presence of sharp corners in single 
wires and the interaction between circular wires result in absorption and energy-loss peaks 
to be broadened by the existence of multipolar surface resonances. 

J. M. P. acknowledges partial support by the Basque Unibertsitate eta Ikerketa Saila 
and the Spanish Ministerio de Education y Cultura. 
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they lose energy Two while undergoing a wave vector change q is proportional to the 
so-called energy-loss function, i.e., the imaginary part of the finite wave-length effective 
inverse longitudinal dielectric function. In typical electron energy-loss experiments the 
energy- loss function can sometimes be approximated by the g — *> limit, and in this 
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and n = 1. 

[33] That the n depolarization factor equals / and is, therefore, independent of the geometry 
of the two-dimensional structure is only valid, within the quasistatic approximation, for 
small values of the length d characterizing our two-dimensional structure. In the case 
of metallic wires of the order of 1/zm radius the self- inductance of the wire structure 
results in the no depolarization factor to decrease, for a given filling fraction, as ~ 1/d 2 
(see J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, Phys. Rev. Lett. 25, 
4773 (1996)). 

[34] For three-dimensional composite structures, D = 3, a combination of Eqs. (2) and (5) 
with Eq. (7) results, if all multipolar modes are neglected, in the dielectric function 
originally derived by Maxwell-Garnett [15] for a system of spherical particles. 



[35] Although the results presented in Figs. 2 and 3 are all well converged for x < 2, the 
results presented for x = 6 (corresponding to a sampling mesh of 120 x 1 x 120) might 
slightly differ at u ~ ui p / y/2 from well-converged results. 
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FIGURES 

FIG. 1. Periodic system of metallic wires with square cross sections of side d, arranged in a 
square array with lattice constant a. The wires are infinitely long in the y direction. 

FIG. 2. Imaginary part of the effective long- wavelength dielectric function of the periodic sys- 
tem described in Fig. 1, for p polarized electromagnetic excitations and various values of the 
volume filling fraction: 25.0% (dotted line), 44.4% (dashed line), and 75.6% (solid line). The 
results obtained for a volume filling fraction of 2.3% are represented in the inset. 

FIG. 3. Same as in Fig. 2, for the effective energy-loss function. The bulk energy-loss function 
(x = 1) is represented by a thick solid line. 

FIG. 4. Dipolar mode positions (depolarization factors), as a function of the volume filling 
fraction, for wires with square cross sections. Solid lines: our numerical results for mi (the line 
below 0.5) and n\ (the line above 0.5). Dotted lines represent the MG depolarization factors given 
by Eqs. (8) and (9). Dipolar mode strengths are represented in the inset: Solid and dashed lines 
represent the coefficients B\ and C\ entering Eqs. (2) and (5), respectively. The MG dipolar mode 
strengths, B\ and C±, are both equal to unity. 
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Figure 1 
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Figure 4 




